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Abstract 


Further  Simplified  Characterizations  of  Linear 
Complementarity  Problems  Solvable  as  Linear  Programs 

Faiz  A.  Al-Khayyal 
Georgia  Institute  of  Technology 


Mangasarian  developed  necessary  and  sufficient  conditions  under 
which  a  linear  complementarity  problem  can  be  solved  as  a  linear  program. 
These  conditions  include  strict  inequality  systems  which  are  difficult  to 
apply.  Using  a  bilinear  programming  approach  and  column  scaling,  we  are 
able  to  derive  further  simplified  conditions  which  replace  the  strict 
inequality  systems  with  equations.  Our  simplified  conditions  yield  new, 

checkable  special  cases  of  linear  complementarity  problems  that  are 

/ 

solvable  as  linear  programs.  ^ - - 

Key  Words:  Linear  Complementarity  Problem,  Linear  Programming,  Bilinear 
Programming 


Abbreviated  Title:  Conditions  for  Solving  LCP  as  LP 


This  work  was  supported  in  part  by  the  Office  of  Naval  Research  under 
Contract  No.  N00014-86-K-0173  and  by  general  research  development  funds 
provided  by  the  Georgia  Institute  of  Technology.  Reproduction  is 
permitted  in  whole  or  in  part  for  any  purpose  of  the  U.S.  Government. 


Further  Simplified  Characterizations  of  Linear 
Complementarity  Problems  Solvable  as  Linear  Programs 

Faiz  A.  Al-Khayyal 
Georgia  Institute  of  Technology 


1.  Introduction 

In  a  series  of  papers  [3, 4, 5,6],  Mangasarian  investigated  necessary 
and  sufficient  conditions  under  which  the  linear  complementarity  problem 
(LCP)  of  finding  a  real  n-vector  x  such  that 

Mx  +  q  >  0,  x  >  0,  x^(Mx  +  q)  =  0  (1) 

is  solvable  (has  a  nonempty  solution  set),  where  M  and  q  are  given  real 
nxn  and  nxl  matrices,  respectively.  For  specified  nonnegative  vectors  r 
and  s  that  appear  in  these  conditions,  the  linear  program  (LP) 

T  T 

minimize  (r  +s  M)  x 
x 

subject  to  Mx  +  q  >  0,  x  >  0  (2) 

has  a  nonempty  optimal  solution  set  which  is  a  subset  of  the  LCP's 
solution  set  (see  [5]  for  details).  Thus,  such  conditions  characterize 
LCP's  that  can  be  solved  as  LP's.  In  a  follow-up  paper  [6],  simpler 
characterizations  were  derived. 

Aside  from  the  obvious  value  of  being  able  to  solve  some  LCP’s 
(which  are  NP-complete)  as  LP's  (which  are  polynomially  solved),  such 
conditions  may  be  useful  in  developing  constructive  characterizations  of 
matrix  classes  that  possess  complementary  solutions  (see,  e.g. ,  [1]), 


*r, 

to 
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especially  subclasses  that  are  not  independent  of  the  vector  q.  In  this 


note,  we  use  a  bilinear  programming  approach  to  give  a  new,  more  direct 


derivation  of  conditions  that  further  simplify  the  characterizations  in 


[6].  We  use  our  conditions  to  derive  new  special  cases  under  which  an 


LCP  is  solvable  as  an  LP.  In  particular,  unlike  most  of  the  cases  in 


[5],  our  simplified  characterizations  allow  for  new  cases  that 


incorporate  the  vector  q. 


The  following  notation  is  used  in  this  paper.  The  space  of  mxn  real 


matrices  is  denoted  by  RmXn.  We  write  Rm**  =  Rra  and  use  a  superscript  T 


to  transpose  all  matrices.  The  vector  in  R  having  ones  for  all  of  its 


elements  is  denoted  by  e  and  IeR  is  the  usual  identity  matrix.  A 


Z-matrix  is  a  square  matrix  with  nonpositive  off-diagonal  elements  and  Z 


is  used  to  denote  the  class  of  such  matrices.  The  class  KcZ  consists  of 


nonsingular  Z-matrices  whose  inverses  are  nonnegative. 


2.  Previous  Work 


The  results  of  Mangasarian  [6]  are  summarized  In  this  section. 


Interested  readers  should  consult  the  source  for  full  proofs. 


The  dual  program  (DP)  of  the  linear  program  (2)  is  defined  as 


T 

maximize  -q  y 

y 


T  T 

subject  to  -H  y  +  r  +  M  s  )  0,  y>0. 


This  problem  plays  an  important  role  in  proving  the  next  three  theorems 


which  characterize  the  conditions  under  which  an  LCP  can  be  solved  as  an 
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Theorem  2.1  The  LCP  has  a  solution  if  and  only  if  the  LP  is  solvable  for 
some  r,s  e  Rn,  r,s  >  0,  such  that  an  associated  dual  optimal 
variable  y  satisfies 

(I-M)Ty  +  r  +  MTs  >  0. 

Furthermore,  each  solution  of  the  LP  solves  the  LCP. 

Theorem  2.2  If  there  exist  c,r,s  e  Rn  and  Z^,Z2  e  Rn*n  such  that 

MZ^  =  +  qc"^,  rTZ^  +  s"*^  >  0,  c,r,s  >  0,  Z^,Z2  e  Z, 

(4) 

then  the  feasible  LCP  has  a  solution  which  can  be  obtained 
by  solving  the  LP.  Conversely,  if  the  LCP  has  at  least  one 
solution  x  which  is  a  vertex  of  {x:  Mx  +  q  >  0,  x  >  0}  which 
is  also  nondegenerate,  that  is,  x+Mx+q>0,  then 
conditions  (4)  are  satisfied  and  each  solution  of  the  LP 
solves  the  LCP. 

Theorem  2.3  The  LCP  has  a  solution  if  and  only  if  the  LP  is  solvable  for 
some  r,s  e  R°  which  must  satisfy  the  following  conditions: 

(a)  MZX  =  z2  +  qcC 

(b)  r1^  +  sTZ2  >  0 

(c)  rTZ1  '+  sTZ2  +  cT  >  0  (5) 


for  some  vector  c  e  Rn  and  some  matrices  Z^,Z ^  e  Rn*n. 
Furthermore,  each  solution  of  the  LP  solves  the  LCP. 

By  examples  Mangarasian  [6]  shows  that  the  nondegenerate  vertex 
assumption  in  Theorem  2  cannot  be  relaxed.  Theorem  1  strengthens  the 
sufficient  conditions  of  [3,  Lemma  1]  to  be  necessary.  Theorem  2 
simplifies  the  sufficient  conditions  of  [4,  Theorem  1]  and  strengthens 
them  to  be  also  necessary  under  the  nondegeneracy  assumption.  Theorem  3 
simplifies  the  sufficient  conditions  of  [5,  Theorem  2]  which  were  shown 
to  be  also  necessary  by  Pang  [7]. 

3.  Some  Simplified  Characterizations 

The  theorems  in  the  preceding  section  give  conditions  under  which 
appropriate  nonnegative  vectors  r  and  s  exist  such  that  the  solution  of 
the  linear  program  (2)  will  also  solve  the  LCP.  If  r  and  s  are  allowed 
to  vary  in  (2),  the  natural  question  to  ask  is  under  what  conditions  will 
minimizing  over  these  variables  along  with  x  produce  a  solution  to  the 
LCP.  Of  course,  in  this  situation  the  objective  function  of  (2)  is 
bilinear.  The  answer  to  the  question  turns  out  to  be  surprisingly 
simple. 

Consider  the  following  bilinear  program  (BLP) 


T  T  T 

minimize  (r  +s  M)x  +  q  s 
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subject  to  Mx+q>0,  x>0 

r  +  s  =  e,  r  >  0,  s  >0. 


Notice  that  for  fixed  (r,s)  the  bilinear  program  (6)  reduces  to  the 
linear  program  (2)  (ignoring  the  feasibility  conditions  on  r  and  s 
above.)  In  addition,  we  now  also  have  a  linear  program  in  (r,s)  for  fixed 


Lemma  3.1  The  LCP  has  a  solution  if  and  only  if  the  BLP  has  a  vanishing 
optimal  objective  value. 

Proof:  Rewrite  (6)  as  the  bilinear  program 


T  T 

minimize  (e-s)  x  +  s  (Mx  +  q) 
(x,s) 


subject  toMx+q>0,  x>0  (7) 

0  <  s  <  e. 

The  proof  follows  from  the  observation  that  the  objective 
function  is  nonnegative  over  the  feasible  region.  (See 
Appendix  for  details.) 


For  fixed  s,  the  bilinear  program  (7)  becomes  the  linear  program 


T  T  T 

minimize  (e-s  +  M  s)  x  +  q  s 


subject  to  Mx  +  q  >0,  x  >  0. 


The  dual  of  (8)  is  the  linear  program 
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maximize  -q  (y-s) 

y 

T 

subject  to  M  (y-s)  <  e-s,  y  >  0. 


Under  the  transformation  y  =  u  +  s,  this  is  equivalent  to  the  problem 


maximize  -q  u 
u 

T 

subject  to  M  u  <  e-s,  u  >  -  s. 


(9) 


When  the  roles  of  (2)  and  (3)  in  the  results  of  Mangasarian  [6]  are 
assumed  by  (8)  and  (9)  above,  the  three  theorems  of  the  preceding  section 
can  be  further  simplified.  By  column  scaling  we  are  also  able  to  replace 
the  strict  inequalities  in  (A)  and  (5c)  by  equations.  Proofs  for  the 
next  three  theorems  are  provided  in  the  Appendix. 

Theorem  3.2  The  LCP  has  a  solution  if  and  only  if  the  linear  program  (8) 
is  solvable  for  some  s  e  Rn,  0  <  s  <  e,  such  that  an 
associated  dual  optimal  variable  u  satisfies 


(M-I)  u  <  e. 


(10) 


Furthermore,  each  solution  of  the  linear  program  (8) 
solves  the  LCP. 


nxn 


Theorem  3.3  If  there  exist  c,s  e  R  and  Z^, e  R  such  that 
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(a)  MZ  «  Z2  +  qc 

(b)  (e-s)TZ1  +  sTZ2  =  eT 

c  >  0,  0  <  s  <  e,  Z1,Z2  e  Z, 


(ID 


then  the  feasible  LCP  has  a  solution  which  can  be  obtained 
by  solving  the  linear  program  (8).  Conversely,  if  the  LCP 
has  at  least  one  solution  x  which  is  a  vertex  of  (x:  Mx  +  q 

>  0,  x  >  0}  which  is  also  nondegenerate,  that  is,  x  +  Mx  +  q 

>  0,  then  the  conditions  (11)  are  satisfied  and  each 
solution  of  the  linear  program  (8)  solves  the  LCP. 

Theorem  3.4  The  LCP  has  a  solution  if  and  only  if  the  linear  program  (8) 
is  solvable  for  some  s  e  Rn  which  must  satisfy  the  following 
conditions: 


(a)  HZl  =  Z2  +  qc 

T  T  T  T 

(b)  (e-s)  Zx  +  s  Z2  +  c  =  e 

0<c<e,  0<s<e,  e  z* 


(12) 


,nxn 


for  some  vector  c  e  R  and  some  matrices  Z^,Z2  e  R' 
Furthermore,  each  solution  of  the  linear  program  (8)  solves 
the  LCP. 
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In  addition  to  the  cases  in  [2]  and  [5],  new  special  cases  can  be 

T 

derived  from  the  conditions  (11)  and  (12).  Let  p  =  e+(M-I)  s  denote  the 
objective  cost  vector  and  let  S  =  {x:  Mx+q  >  0,  x>0)  denote  the  feasible 
set  of  the  linear  program  (8).  The  following  theorem  provides  new  cases 
such  that  the  LCP  is  solvable  as  a  linear  program.  In  particular,  the 
simplified  conditions  (11)  and  (12)  allow  for  the  easy  incorporation  of 
the  vector  q  into  our  cases. 


Theorem  3.5 


Let  S*0,  and  let  (M,q)  satisfy  any  of  the  conditions  below. 
Then  the  linear  complementarity  problem  (1)  has  a  nonempty 
solution  set,  at  least  one  element  of  which  can  be  obtained 
by  solving  the  linear  program  (8)  with  the  p  indicated. 


(a) 

M-qcteZ,  0<c<e,  (I-M+qcT)T 

s  =  c 

:,  0 

<s<e,  p  = 

e+(M-I)^s 

(b) 

M-qcTeZ,  c>0,  (I-M+qcT)Ts 

=  o. 

0<s 

<e,  p  =  e+(M-I)Ts. 

(c) 

-M-qcTcZ,  0<c<e,  (I-M-qcT) 

T 

s  = 

2e- 

c,  0<s<e, 

p  =  e+(M-I)^s. 

(d) 

-M-qcTeZ,  c>0,  (I-M-qc^)Ts 

=  2e 

,  o 

<s<e,  p  = 

e+(M-I)Ts 

(e) 

T 

=  Z2+qc  j  Z^eZ, 

T 

Zl6  = 

e» 

c!>0,  p  = 

e. 

(f) 

MZX  =  Z2+qcT,  Z  eZ,  Z2eZ, 

T 

Zl6  " 

e- 

c,  0<c<e, 

p  =  e. 

(g) 

MZ^  =  Z2+qc^,  Z^eZ,  Z2eZ, 

T 

Z2e  = 

e. 

c>0,  p  - 

T 

Me. 

(h) 

MZX  =  Z2+qcT,  Zl£Z,  Z2eZ, 

T 

Z2e  = 

e- 

c,  0<c<e, 

T 

p  =  Me. 

(i) 

(M-I)~1eZ,  0<eT(M-I)_1<eT, 

P  = 

(M- 

I)Te. 

(j) 

(I-M)-1eZ,  eT<eT(I-M)-1<2e 

T 

»  P 

=  ( 

I-M)Te. 

(a) 

In  Theorem  3.4,  set  Z^  =  I 

and 

Z2 

=  M-qcT. 

Proof: 


t 

t 


(b) 

In 

Theorem 

3.3, 

set 

=  I  and  Z 2  =  M-qc  . 

(c) 

In 

Theorem 

3.4, 

set 

Z^  =  -I  and  Z2  = 

T 

-M-qc  . 

(d) 

In 

Theorem 

3.3, 

set 

Z^  =  -I  and  Z2  = 

T 

-M-qc  . 

(e) 

In 

Theorem 

3.3, 

set 

s  =  0. 

(f) 

In 

Theorem 

3.4, 

set 

s  =  0. 

(g) 

In 

Theorem 

3.3, 

set 

s  =  e. 

(h) 

In 

Theorem 

3.4, 

set 

s  =  e. 

(i) 

In 

Theorem 

3.3, 

set 

Z1  =  (M— I)-1,  Z2 

=  I+(M-I)_1, 

c  = 

s  0,  and 

T 

s  = 

eT[l-(M-I)-1]. 

(j) 

In 

Theorem 

3.3, 

set 

Z1  =  (I-M)"1,  Z2 

=  (I-M)'1-!, 

c  =  0,  and 

T 

s  = 

eT[(I-M)-1-l]. 

Remark.  A  given  LCP  may  fall  in  more  than  one  case  above,  each  of  which 
could  lead  to  a  different  complementary  solution.  Cases  (a)  and  (b)  are 
not  subsets  of  each  other  since  we  can  construct  problems  that  fall  in 
either  and  in  both.  The  same  holds  for  cases  (c)  and  (d),  (e)  and  (f), 
and  (g)  and  (h).  In  part  (b),  we  can  always  take  s=0,  but  other  values 
may  work  in  some  instances  so  alternative  solutions  could  be  found. 
Excluding  (e)  through  (h) ,  the  remaining  cases  are  easily  checkable. 

Other  simplified  conditions  can  be  derived  from  the  above  theorem  as 
illustrated  by  the  following  corollaries. 

Corollary  1  Let  S*0,  and  let  M  satisfy  the  conditions  below.  Then  the 


LCP  (1)  has  a  solution  that  can  be  obtained  by  solving  the 
LP  (8)  with  the  p  indicated. 


Proof : 


(a)  In  parts  (a)  and  (b)  of  Theorem  3.5,  set  s  =  c  =  0 


Also  can  be  obtained  from  parts  (e)  and  (f)  by  setting 
c  =*  0  and  *  I. 

(b)  Assume  M  is  nonsingular  and  set  Z^  *  I  and  c  =  0  in 
either  parts  (g)  or  (h)  of  Theorem  3.5. 

(c)  Subsumed  by  part  (i)  of  Theorem  3.5. 

(d)  Subsumed  by  part  (j)  of  Theorem  3.5. 

(e)  In  parts  (c)  or  (d)  of  Theorem  3.5,  set  c  =  0. 

(f)  In  parts  (g)  or  (h)  of  Theorem  3.5,  set  Z^  =  -I  and 

c  =  0. 

Corollary  2  Let  S*0,  and  let  (M,q)  satisfy  the  conditions  below.  Then 
the  LCP  (1)  has  a  solution  that  can  be  obtained  by  solving 
the  LP  (8)  with  the  p  indicated. 

(a)  M  =  -qeT,  p  “  M^e. 

T  T 

(b)  M  =  qs  ,  0<s<e,  p  =  e+(M-I)  s. 

T 

(c)  M  =  qc  ,  c>0,  p  =  e. 

T  T 

(d)  M  *  -I-qc  ,  c>0,  p  *»  M  e. 

T 

(e)  M  *  -I+qc  ,  c>0,  p  ■  e. 

T  T 

(f)  M  ■  I+qc  ,  c>0,  p  ■  e+(M-I)  s,  where  0<s<e. 

T 

(g)  M-qc  eZ,  c>0,  p  -  e. 

(h)  M  *(I+qc^)eZ,  c>0,  p  ■  M^e. 


Proof : 


In  Theorem  3.5,  set: 


(a)  Zj  -  -I,  Z2  -  0,  and  c=e  in  part  (h). 

(b)  c  =  s  in  part  (a). 

(c)  Z1  =  I  and  Z2  =  0  in  part  (e). 

(d)  Z^  =  -I  and  Z2  =  I  in  part  (g). 

(e)  Z^  =  I  and  Z2  -  -I  in  part  (e). 

(f)  I-M+qcT  =  0  in  part  (b). 

(g)  s  =  0  in  part  (b). 

(h)  Z1  =  M_1(I+qcT)  and  Z2  =  I  in  part  (g). 


Remark.  Cases  (a)  and  (b)  of  Corollary  1  are  well-known.  They  are  also 
special  cases  of  (g)  and  (h),  respectively,  in  Corollary  2.  Notice  that 
M-qc  eZ  for  some  c  >  0  if  and  only  if  for  each  j  *  l,...,n,  either  < 
0  for  all  i*j  or  <  c^q  for  all  i*j  and  some  c^X),  where  is  the 
(i.j)-element  of  M  and  c^  is  the  jth  component  of  c.  In  other  words, 
excluding  the  diagonal  elements,  if  every  column  of  M  is  either 
nonpositive  or  does  not  exceed  a  positive  multiple  of  q,  then  a  solution 
of  the  associated  LCP  can  be  obtained  by  solving  a  linear  program  with 


objective  cost  vector  e.  In  particular,  a  problem  does  not  satisfy  this 
condition  if  an  off-diagonal  element  of  M  is  positive  and  the  element  of 
q  in  the  same  row  is  nonpositive.  On  the  other  hand,  we  can  easily 
characterize  all  vectors  q  that  do  satisfy  this  condition. 
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Appendix 

Proof  of  Lemma  3.1 

Suppose  x  solves  the  LCP  and  let  s  denote  the  n-vector  whose  ith 
component  is  defined  as 


fl  if  xt  >  0 

0  if  ^  =  0 


(13) 


Clearly,  (x,s)  solves  (7)  with  a  vanishing  optimal  objective  value. 

Now  suppose  that  (x,s)  solves  (7)  and  the  objective  function 
vanishes  at  this  point.  Then 


(l-si)xi  =  0  and  s^Mx  +  q^  =  0 

-T  - 

for  all  i  =  l,2,...,n,  which  implies  that  x  (Mx  +  q)  =  0. 

Proof  of  Theorem  3.2 

Suppose  x  solves  the  LCP  and  let  s  be  specified  by  (13).  By  Lemma 

-  t- 

3.1,  x  solves  (8)  for  s  =  s  and  q  u  =  0  for  all  optimal  solutions  u  of 
(9).  In  particular,  u  =  0  solves  (9)  and  satisfies  (10). 

For  some  0  <  s  <  e,  suppose  that  x  solves  (8)  and  its  complementary 
dual  solution  u,  which  solves  (9),  satisfies  (10).  By  complementary 
slackness,  we  have 

x^(e-s  -  M^u)  “  0  and  (u+s)^(Mx  +  q)  ■  0. 


If  x  does  not  solve  the  LCP,  then  there  is  an  i  such  that  x^(Mx  +  q)^  >  0 
so  that 


(e-s  +  M^u) 


0  and  (u  +  s). 


0. 


T-  _ 

Adding  gives  [e-(I-M)  u]  ■  0  which  contradicts  the  assumption  that  u 

satisfies  (10).  Hence,  x  solves  the  LCP.  To  show  the  second  part  of  the 

A 

theorem,  let  (x,  u,  s)  be  as  above  and  let  x  t  x  also  solve  (8)  for 
s  =  s.  By  the  strong  duality  theorem  we  have 


0 


(e-s+M^s) ^x  +  q^s  +  q^u  +  (M^u)^x  -  (MTu)^x 
(e-s  -  M^u)^x  +  (u+s)^(Mx  +  q). 


Since  each  term  is  nonnegative,  we  must  have 

*rp  w  rp  A 

x  (e  -  s  -  Mu)  =  0  and  (u+s)  (Mx  +  q)  ■  0. 

A 

Because  u  satisfies  (10),  we  conclude  by  our  argument  above  that  x  solves 

a 

the  LCP.  Notice  that  the  complementary  dual  solution  of  x  need  not 
satisfy  (10). 

Proof  of  Theorem  3.3 

Assume  the  LCP  is  feasible,  let  (s,  c,  Z^ ,  Z solve  (11)  and  let  x 
solve  (8)  for  s  -  s.  Then  x  solves  (2)  for  (r,  s),  where  r  *  e  -  s,  and 
(r,  s,  c,  Z j,  zp  solves  (4).  Hence,  by  Theorem  2.2,  the  LCP  has 


a 
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solution  that  can  be  obtained  by  solving  (2),  which  is  equivalent  to  (9) 
when  r  ■  e  -  s. 

Now  suppose  x  is  a  nondegenerate  extreme  point  solution  of  the  LCP 
and  let  s  be  specified  by  (13).  By  Lemma  3.1,  x  solves  (8)  for  s  -  s 
with  vanishing  optimal  objective  value.  Hence,  u  *  0  is  the  unique 
solution  of  (9).  Since  u  satisfies  (10),  by  Theorem  3.2,  each  solution 
of  (8)  solves  the  LCP.  Moreover,  conditions  (11)  are  satisfied  since, 
following  the  proof  of  Theorem  2.2  in  [6],  they  are  equivalent  to 


\ 

T 

M  u  <  e  -  s 

u  >  -  s 
T 

q  u  <  l 


->  (M  -  I)Tu  <  e 


J 

for  some  o  <  s  <  e,  which  is  clearly  satisfied  by  (s,  u)  above. 

Proof  of  Theorem  3.4 

Suppose  x  solves  the  LCP  and  let  s  be  specified  by  (13).  By  Lemma 
3.1,  x  solves  (8)  for  s  -  s.  Let 


!*  Z1 


-  T  — 
-  xe  and  Z. 


-  (Mx  +  q)e 


Then  (s,  c,  Z^,  Z^)  satisfies  (12). 

Now  let  (s,  c,  Z  t  2^)  satisfy  (12)  and  let  x  solve  (8)  for  s  *  s. 
With  r  ■  e  -  s,  it  is  easy  to  verify  that  (r,  s,  c,  2^  satisfies  (5) 
and  x  solves  (2)  for  (r,  s).  Hence,  by  Theorem  2.3,  the  LCP  has  a 
solution  and  each  solution  of  (2),  which  is  equivalent  to  (8)  when 


r  ■  e  -  s,  solves  the  LCP 
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